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Abstract 

Paired sense and antisense (S/AS) genes located in as represent a structural feature common to the genomes of both 
prokaryotes and eukaryotes, and produce partially complementary transcripts. We used published genome and 
transcriptome sequence data and found that over 20% of genes (645 pairs) in the budding yeast Saccharomyces cerevisiae 
genome are arranged in convergent pairs with overlapping 3'-UTRs. Using published microarray transcriptome data from 
the standard laboratory strain of S. cerevisiae, our analysis revealed that expression levels of convergent pairs are 
significantly negatively correlated across a broad range of environments. This implies an important role for convergent 
genes in the regulation of gene expression, which may compensate for the absence of RNA-dependent mechanisms such as 
micro RNAs in budding yeast. We selected four representative convergent gene pairs and used expression assays in wild 
type yeast and its genetically modified strains to explore the underlying patterns of gene expression. Results showed that 
convergent genes are reciprocally regulated in yeast populations and in single cells, whereby an increase in expression of 
one gene produces a decrease in the expression of the other, and vice-versa. Time course analysis of the cell cycle illustrated 
the functional significance of this relationship for the three pairs with relevant functional roles. Furthermore, a series of 
genetic modifications revealed that the 3'-UTR sequence plays an essential causal role in mediating transcriptional 
interference, which requires neither the sequence of the open reading frame nor the translation of fully functional proteins. 
More importantly, transcriptional interference persisted even when one of the convergent genes was expressed ectopically 
(in trans) and therefore does not depend on the cis arrangement of convergent genes; we conclude that the mechanism of 
transcriptional interference cannot be explained by the transcriptional collision model, which postulates a clash between 
simultaneous transcriptional processes occurring on opposite DNA strands. 
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Introduction 

Sense and anti-sense transcripts (S/AS) are simply pairs of 
RNAs (either protein-coding or non-protein-coding) containing 
sequences that are at least partially complementary to each other. 
S/AS gene pairs can be transcribed in cis from opposing DNA 
strands at the same genomic locus [1]. In recent years, genome 
sequencing projects have revealed the frequent presence of anti- 
sense transcripts in cells, even in one of the smallest self-replicating 
organisms. Mycoplasma pneumoniae [2]. In particular, S/AS pairs 
located in cis represent a common structural feature in the genomes 
of both prokaryotes and eukaryotes, including mammals (human, 
mouse, rat, cow), birds (chicken), lower vertebrates (zebrafish), 
invertebrates [Camorhabditis elegans, Drosophila melanogaster), plants 
{Arabidopsis thaliana, rice), and yeast [3-6] . According to their tran- 
scriptional orientation and extent of sequence overlap, S/AS genes 
can be classified into three major groups: (1) convergent gene parrs, 
overlapping at their 3' ends; (2) divergent gene pairs, overlapping at 
their 5' ends; and (3) consistent gene pairs, overlapping and tran- 
scribed in the same direction. The genomic distribution of these 



different types of gene pair is species-specific. For example, con- 
vergent gene pairs are prevalent in Drosophila and C. elegans, but are 
rare in human and mouse genomes [3] . 

The structural organization of S/AS gene pairs confers a 
significant mechanism for regulating gene expression levels. For 
example, global analysis of the mammalian transcriptome showed 
that a large proportion of the genome can produce transcripts 
from both strands and revealed anti-regulation of S/AS pairs [7]. 
In particular, experimental perturbation of an antisense RNA was 
shown to alter the expression of the sense mRNA in cis. In both 
budding and fission yeasts, transcriptional interference has been 
observed between sense and antisense transcripts of several 
convergent genes arranged in cis or trans [8-12]. This involves 
one transcriptional process exerting a direct negative impact on a 
second transcriptional process. Such interference plays important 
functional roles, for example in the entry of yeast cells into meiosis 
by regulating the expression of the IME4 gene [9] . 

Antisense transcripts are mature RNA species (polyadenylated 
at their 3' ends) and contribute a mechanism of transcription 
interference that is distinct from the RNAi-mediated regulation of 
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Author Summary 

In the compact genonne of the budding yeast Saccharo- 
myces cerevisiae, genes are frequently organized into 
convergent pairs that are transcribed from opposing 
DNA strands in opposite directions and have overlapping 
3'-UTRs. Here we explore the negative correlation in 
expression levels between convergent genes using a set of 
645 convergent pairs in the yeast genome, Identified from 
published genomic and transcrlptomic sequence data and 
accounting for ~20% of total yeast genes. Analysis of 
published microarray experiments confirmed that the 
negative correlation In expression between convergent 
genes occurs across a broad range of growth conditions. 
This Implies that such transcriptional interference is an 
Important means of regulating gene expression in yeast, 
especially in the absence of other eukaryotic RNA- 
dependent mechanisms such as micro RNAs. We focused 
on profiling the expression of four pairs of convergent 
genes In wild type yeast and Its genetically modified 
strains, to explore the causes and mechanisms of 
transcriptional interference. We demonstrate that the 3'- 
UTR sequence alone plays the essential and causal role in 
Interference between convergent genes. Intrlguingly, 
transcriptional interference occurs even when one of the 
convergent genes Is expressed from elsewhere in the 
genome (In trans), raising new questions about how 
transcriptional Interference operates. 



gene expre.ssion that Ls present in most eukaryotes, but absent in 
Saccharomyces cerevisiae. For example, antisense transcripts involved 
in RNAi are often encoded elsewhere in the genome from the target 
gene, and require further processing into shorter, functional 
sequences by the RNAi machinery. The most obvious explanation 
for transcriptional interference between convergent S/AS gene pairs 
is given by the coUision model [13]; in this model, RNA synthesis 
from one DNA strand clashes with transcription from the other 
strand, and so active antisense transcription would suppress sense 
RNA transcription [14,15]. Though this model is supported by 
atomic force microscopy data in E. coli [16], its role in transcriptional 
interference in budding yeast has not been thoroughly assessed. 

In this work, we have exploited the availability of genome-wide 
transcriptome data from the budding yeast S. cerevisiae to explore 
the significance of transcriptional interference between convergent 
gene pairs on a global scale. Using a comprehensive set of S/AS 
pairs with overlapping 3'-UTRs, we have demonstrated that tran- 
scriptional interference is common to more than 600 gene pairs, 
which account for ~20% ORFs in the S. cerevisiae genome, across a 
broad range of growth conditions. For a more detailed under- 
standing of the underlying mechanisms, we focus on four such 
representative gene pairs and show that transcriptional interfer- 
ence is dependent only on the 3'-UTR sequence. Moreover, its 
occurrence is not restricted to the arrangement of the gene pairs in 
cis, but can also occur when the partner genes are re-located apart 
(in trans). We have illustrated the functional importance of this 
mode of gene regulation during the yeast cell cycle and its role in 
the phenotypic response of yeast cells to environmental stress. 

Results 

Genome-wide transcriptional interference of convergent 
genes in the yeast genome 

Based on genome-wide transcriptome [17] and genome sequence 
data, we identified 645 convergent gene (ORF) pairs with overlapping 



3'-UTRs, accounting for ~20% of total genes in the yeast genome 
(see materials and methods. Table SI). We compared the con- 
vergent ORF pairs predicted from the mRNA-Seq data with those 
from both nascent RNA sequencing (NET-Seq) data [18] and 
strand-specific RNA sequencing (ssRNA-Seq) data [19]. We found 
that among the 645 convergent pairs, 531 (82.3%) and 329 (51.0%) 
were detected as convergent parrs with overlapping 3'-UTRs from 
the NET-Seq and ssRNA-Seq datasets respectively. For 15.5% and 
33.0% of the 645 convergent pairs, at least one of the two genes in 
the pair was unexpressed, in the NET-Seq and ssRNA-Seq datasets 
respectively (Figure SI). Additionally, there are only 12M sequenc- 
ing reads in the ssRNA-Seq dataset, whilst the mRNA-Seq and 
NET-Seq datasets contain 29M and 69M reads respectively. This 
may, at least pardy, explain the lower proportion of the convergent 
pairs confirmed in the ssRNA-Seq data compared to the NET-Seq 
data. 

Based on the mRNA-Seq and genomic sequencing data, we 
observed a highly significant negative correlation between 3' -UTR 
length and gene expression level of the corresponding ORF 
(»■= —0.089, P<\0 *), suggesting an important role for the 3'- 
UTR in the regulation of transcription. However, there was no 
clear relationship between 5'-UTR length and gene expression level. 
To assess the occurrence of transcriptional interference between 
convergent partner genes with overlapping 3'-UTRs, we extracted 
niRNA expression values for these 645 gene parrs from Affymetrix 
microarray experiments based on the test strain BY4741 under 
seven environmental stress conditions (Table S3). The expression 
of convergent ORFs showed a consistent and highly significant 
negative correlation (-0.159<r<-0.124; 10"^<P<0.002) across 
all seven conditions (Figure S2). The significant negative 
correlation was also observed in both the NET-Seq (r= —0.09, 
P<0.05) and ssRNA-Seq data (r= -0.121, P<0.05). We infer that 
transcriptional interference between convergent gene pairs with 
overlapping 3'-UTRs is a widespread phenomenon in the yeast 
genome. 

Transcriptional interference between convergent genes 
located in cis and in trans 

To further understand the transcriptional relationship between 
genes with overlapping 3'-UTRs, we selected four convergent 
gene pairs from the complete set of 645 pairs for more detailed 
study (Table 1). The 3'-UTR boundaries and overlapping regions 
were confirmed by 3 '-RACE sequencing. In addition, these four 
pairs were also confirmed as convergent genes with overlapping 
3'-UTRs in the NET-Seq and ssRNA-Seq datasets (Figure S3). 
We profiled gene expression of the convergent pairs in the wild 
type strain YLIC^^, a laboratory haploid strain which we have 
described previously [20], and in a series of seven derived genet- 
ically engineered strains (Table 2), as illustrated in Figure lA. 
When expression of upstream genes [SHMl, AXL2, APTl and 
ADEl) was inhibited (group I), the expression of the downstream 
partner genes [YPTIO, REV7, UJVGl and KIN3) was up-regulated 
in comparison to the wild type control. Correspondingly, when the 
downstream genes were over-expressed, expression of the 
upstream partner genes was clearly repressed (group II). In the 
extreme situation where expression of the downstream genes was 
completely silenced (group III), expression of the upstream genes 
increased to a variable extent compared to the wild type level. 
Convergent overlapping gene pairs therefore exhibit a pattern 
whereby a change in the expression of either partner gene leads to 
the expression of the other gene changing in the opposite direc- 
tion, a pattern we refer to as 'anti-regulation' hereafter. 

To explore the relationship between transcriptional interference 
and the sequence of the ORF, and thus of the protein encoded, we 
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Table 1. Four convergent pairs of S. cerevislae ORFs with overlapping 3'-UTRs. 





ORF pair 


3'-UTR length 


Length of sequence overlap 


Biological process 



AXL2 


127 bp 


164 bp 


Axial cellular bud site selection In cell-cycle 


REV7 


178 bp 




Error-free/error-prone translesion synthesis 


APTl 


211 bp 


157 bp 


Adenine biosynthesis 


UNGl 


180 bp 




DNA repair, regulated by cell-cycle 


ADEl 


98 bp 


206 bp 


Purine nucleotide biosynthetic process 


KIN3 


206 bp 




Protein kinase, mitotic spindle orientation 


SHMl 


53 bp 


29 bp 


Serine family amino acid biosynthesis 


YPTIO 


9 bp 




Golgi organization 



The upstream (downstream) gene is shown in row 1 (2). 
doi:l 0.1 371 /journal.pgen.l 004021 .tOOl 



created frame-shift mutations in the coding sequence of the 
downstream genes. These mutations were designed to generate 
premature termination in translation, which, while not changing 
the normal transcription of the genes, would mean that the 
transcripts generated would no longer be translated into normally 
functioning proteins. However, anti-regulation (transcriptional 
interference) was still observed in the presence of such mutations 
(group IV). Expression of the frame-shifted downstream genes was 
increased relative to that of the wild types, while at the same time, 
expression of the upstream genes was correspondingly decreased. 
This shows that over-expressing the mutated downstream genes 
leads to significantly altered expression of their upstream conver- 
gent partners, and thus excludes the possibility that anti-regulation 
of gene partners in these four convergent pairs is dependent on 
functional interactions between their encoded proteins. 

Having established the transcriptional interference between the 
four pairs of convergent genes with overlapping 3'-UTRs, we then 
created a series of genetically modified strains to test if the 
interference could be maintained in trans. These modifications 
included translocation of either the complete (ORF and 3' -UTR) 
or partial (ORF or 3' -UTR only) sequence of the downstream 
genes to the HO locus and ectopic expression stimulated by the 
constitutive promoter Padh- Ectopic over-expression of the com- 
plete downstream genes (group V) or only their 3'-UTRs (group 
VII) significantly repressed the expression of the upstream genes. 
In contrast, over-expression of the ORF alone did not lead to an 
obvious decrease in expression of the upstream genes (group VI). 
Firstly, this shows that the 3' -UTR (rather than the ORF) plays a 
key role in anti-regulation between the partner genes of convergent 
gene pairs. Secondly, the 3'-UTRs can mediate anti-regulation of 
expression of the partner genes in tram, i.e. when the partner genes 
are located apart, and therefore does not depend on their native 
arrangement in cis. It should be noted that the range of variation in 
relative expression of downstream genes is large in comparison to 
that of the corresponding upstream genes in groups II, IV, V, VI 
and VII. This reflects the fact that we designed the experiments to 
overexpress the downstream genes. 

We compared the transcriptional interference at both mature 
RNA and nascent RNA levels for two of the four convergent 
gene pairs. Figure IB illustrates expression of one gene when 
expression of its convergent partner was altered for two convergent 
gene pairs at nascent RNA (lower panel) and mature mRNA 
(upper panel) levels. It is clear that expression of one of the 
convergent genes markedly responds to altered expression of its 
convergent partner at both nascent and mature RNA levels. This 
excludes the possibility that the anti-regulation of expression 
between convergent gene pairs occurs during the RNA maturation 
process. 



Transcriptional interference between convergent genes 
is reflected in protein abundance 

We next sought to determine whether the observed anti- 
regulation in expression between convergent gene pairs with 
overlapping 3'-UTRs could also be observed at the level of protein 
abundance. We replaced the ORFs of the upstream partners for 
two of the four gene pairs [APT]/ UMGl and ADEl /KIMS) witii die 
reporter gene Flue, and quantified expression of the reporter 
protein in the wild type and modified backgrounds. Over- 
expression of either the downstream genes in situ, or of their 3'- 
UTRs ectopicaUy, caused significant reduction in levels of the 
reporter protein (Figure 2). However, ectopic over-expression of 
only the ORF did not alter levels of the reporter protein. 
Therefore interference between convergent genes with overlap- 
ping 3'-UTRs affects both the expression of gene transcripts and 
protein levels. These effects are dependent on the 3' -UTR 
sequence and can be observed in tran.s as well as in eis. 

To further confirm the causal role of the overlapping 3'-UTRs 
in the transcriptional interference between convergent gene pairs, 
we removed the overlapping 3'-UTRs of two convergent pairs 
[AVTl/UNGl and ADEl /KIMS) and measured their expression 
responses (Figure 3, group I). Inhibition of expression of the 
upstream gene did not result in up-regulation of the corresponding 
downstream partner gene (group II). Similarly, over-expression of 
the downstream gene did not lead to repression of the upstream 
partner gene (group III). These results reveal that the 3'-UTR 
plays a crucial, causal role in the anti-regulation of expression 
between convergent gene pairs with overlapping 3'-UTRs. 

Convergent genes reach peak expression values at 
different times in the cell-cycle 

Three of the four convergent gene pairs involve a biosynthesis 
pathway gene and a cell-cycle dependent gene (Table 1), 
providing an opportunity to investigate the dynamic change in 
expression of the convergent gene pairs over the course of the cell 
cycle. We profiled the expression of each pair over the 100 minute 
time span of the cell cycle and found a consistent pattern whereby 
one ORF reaches its expression peak, while its convergent partner 
is suppressed (Figure 4A-C). The fourth gene pair involves a 
biosynthesis gene [SHMl) and a GTP-binding protein (TPTIO). As 
expected, this pair did not follow the same pattern (Figure 4D). 
This demonstrates the functional importance of transcriptional 
interference between convergent gene pairs in the cell cycle of 
budding yeast. 

We explored the cell-cycle expression pattern of one of the four 
gene pairs with overlapping 3'-UTRs [ADE1/KIN3) when its 
orientation was switched from 'convergent' to 'tandem'. We created 
two genetically modified strains through homologous recombination. 
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Table 2. Summary of genetically modified strains of the standard wild type yeast YLIC. 





Group Modification targeted 


Strains 


Modifications made 


1 Upstream gene inhibited 


S-Yl 


sbml promoter :: Tjef (terminator inserted upstream of SHM/™^) 




X-Rl 


axl2 promoter :: Tjef (terminator inserted upstream of AXL2°^'') 




P-U1 


apt! promoter :: Tjef (terminator inserted upstream of APTl^^'') 




D-Kl 


adel promoter :: T^ef (terminator inserted upstream of ADEl^^'') 


II Downstream gene over- expressed 


S-Y2 


Padh-VPT;0'*' (Padh inserted upstream of YPTKp"'') 




X-R2 


Padh-REV/"' (Padh inserted upstream of REV/""") 




P-U2 


Padh-UNG;"" (Padh inserted upstream of UNG1°"'') 




D-K2 


Padh-KINB"' (Padh inserted upstream of KIN3°"'') 


III Downstream gene knocl<ed-out 


S-Y3 


yptW°'"' ■:. NAT! (YPTW°"'' knocked out) 




X-R3 


rev7°"'' ■:. NATl {REV7°"^ knocked out) 




P-U3 


ungl°'"' :: NATl {UNG1°"'' knocked out) 




D-K3 


kin3°"^ ■:. NATl {KIN3°"'' knocked out) 


IV Downstream gene frame shifted 


S-Y4 


PADH"ypt?0^"' (Frameshift gene yptlO over-expressed by Padh) 




X-R4 


PADH"rev/""' (Frameshift gene rev7 over-expressed by Padh) 




P-U4 


PADH-tJr7g/'""' (Frameshift gene ungl over-expressed by Padh) 




D-K4 


Padh-^'h.^""'* (Frameshift gene I<in3 over-expressed by Padh) 


V Downstream gene (ORF+3'- UTR} ectopically expressed 


S-Y5 


ho ..ypjjoORM'-"™ (YPT10°"'' transferred to HO locus with its 3'-UTR) 




X-R5 


hov.REVf""'*^''"™ (R£l/7°'"' transferred to HO locus with its 3'-UTR) 




P-U5 


ho::UNG;™™' "™ (UNG1°"^ transferred to HO locus with its 3'-UTR) 




D-K5 


ho::KIN3°"''^'''"™ (WNJ™'' transferred to HO locus with its 3'-UTR) 


VI Downstream gene (ORF only) ectopically expressed 


S-Y6 


ho ■YPTIOP"'' {YPTIO""'' transferred to HO locus without its 3'-UTR) 




X-R6 


ho:: REVf""" (REVf""" transferred to HO locus without its 3'-UTR) 




P-U6 


ho::UNGl°"'^ (UNG)™'' transferred to HO locus without its 3'-UTR) 




D-K6 


ho::KIN3°"'' (KIN3°'"' transferred to HO locus without its 3'-UTR) 


VII Downstream gene (3'-UTR only) ectopically expressed 


S-Y7 


ho :: Sh fa/e-l-3'-UTR of yptlO (fusion gene transferred to HO locus) 




X-R7 


ho:: Sh fa/e-t-3'-UTR of rev7 (fusion gene transferred to HO locus) 




P-U7 


ho:: Sh b/e-i-3'-UTR of ungl (fusion gene transferred to HO locus) 




D-K7 


ho:: Sh fa/e4-3'-UTR of I<in3 (fusion gene transferred to HO locus) 


VIII Terminator changed 


P-U8 


Tadh-''IP'"''" & ff.oxrUNGl'" (Changing intergenic region with 
independent terminators) 




D-K8 


"^ADH-^DEl"^ & TAoxr^/W3"'^ (Changing intergenic region with 
independent terminators) 


IX Terminator changed and upstream gene inhibited 


P-U9 


aptl promoter :: Tjef (terminator inserted upstream of APT 1°'^^ in P-U8) 




D-K9 


adel promoter :: Tjef (terminator inserted upstream of ADEl^'^^ in D-K8) 


X Terminator changed and downstream gene over- expressed 


P-UlO 


Padh-UWG;"" (Padh inserted upstream of UNG1°"^ in P-U8) 




D-KIC 


Pf.Dn-K!N3"" (Padh inserted upstream of KIN3°"'' in D-K8) 


Padh is a strong constitutive promoter. 
doi:l 0.1 371 /journal.pgen.l 004021 .1002 



the tandem KIJV3 and the tandem URA3, as illustrated in 
Figure 5A. The modifications were confirmed by sequencing 
(the sequence data is listed in Table S5). Figure 5B illustrates 
the cell cycle expression of the three genes ADEl, KIN3 and 
URA3 in the wild type strain. It is clear that cell-cycle expression 
of the tested convergent gene pair was repeatedly confirmed in 
this independent assay, showing the expression pattern of one- 
rising and the other falling, as shown in Figure 4A. When the 
convergent pair was converted into 'tandem' orientation, the 
ceU-cycle expression patterns of the convergent pair changed in 
comparison to that of the genes in wild type background, but 
the pattern of one rising and the other falling remained 
(Figure 5C). Moreover, we observed the expression pattern of 
ADEl was also markedly altered when its downstream partner 



gene was replaced with a tandemly oriented gene, URA3 
(Figure 5D). These observations strongly support the anti- 
regulation in the cell cycle expression of the convergent pair, 
which is independent of genomic orientation of the convergent 
genes. A natural question arises whether each gene in the con- 
vergent gene pairs tested is expressed in the same cell or instead 
is expressed mainly in different cells in the cell population. To 
answer this question we profiled expression of the convergent 
pair, KIN3 and ADEl, in single cells across different stages of the 
cell cycle (Figure S4). It shows a highly significantly negative 
correlation in expression between the convergent genes (Pear- 
son's correlation coefficient r= —0.38, P<0.01), indicating that 
anti-regulation between convergent genes is also observed within 
individual cells. 
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• SHM1 & YPT10 

• AXL2 & REV7 

• APT1 & UNG1 

• ADE1 & KIN 3 



i 



i 



I 



4 



I 



Wild Type 



III 



IV 



V 



VI 



VII 



B 



• APT1 

• ADE1 



i 



Wild Type 



V 



VI 



2.5 p • UNG1 
• KIN3 

2.0' 

1.5 

1.0 

0.5 

0.0 

-0.5 
1.8 

1.5 

1.2 

0.9 

0.6 

Wild Type 



List of modifications: 

I; Upstream inhibited 
II: Downstream 

overexpressed 
III: Downstream 

knocked-out 
IV: Downstream 

frame-shifted 
V: Downstream 

gene relocated 
VI: Downstream 

only ORF re-located 
VII:Downstream 

only 3'UTR re-located 



Figure 1. RT-PCR expression assay of convergent gene pairs located in cis or in trans. (A) Expression levels of four pairs of convergent 
genes with overlapping 3'-UTRs in wild type yeast, or when the upstream (group I) or the downstream (groups ll-VII) genes were modified by: (I) 
inhibition of gene expression; (II) over-expression; (III) gene knock-out; (IV) frame-shift of the ORF; (V) re-location of the entire ORF and 3'-UTR; (VI) re- 
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location of only the ORF; and (VII) re-location of only the 3'-UTR. The mean and standard deviation of expression level for three replicates are given in 
relative units compared to the wild type genes, which are assigned a value of 1 .0. (B) Expression levels of nascent and mature messenger RNA for one 
of two convergent genes when its convergent partner's expression was altered. 
doi:1 0.1 371 /journal.pgen.1 004021. gOOl 



Phenotypic response to environmental change mediated 
by the 3'-UTRs of convergent genes 

It is well established that gene expression in the budding yeast 
Saccharomyces cerevisiae changes in response to nutrient availability 
[8,21]. For example, ADEl expression is responsive to adenine 
availability in the culture medium; when yeast cells are cultured in 
adenine rich (synthetic complete, SG) medium, the expression of 
ADEl is suppressed, while its expression is stimulated in an 
adenine barren medium (SC-A) [22]. We examined expression of 
the convergent pair ADEl / KIMS in the modified strains D-Kl and 
D-K2 (Table 2), cultured in either adenine rich or barren 
medium. In both strains, ADEl expression was suppressed in the 
SC medium, but enhanced in the SC-A medium. Conversely, 
KIM3 expression was increased in the SC medium and decreased 
in SC-A, even in the presence of the strong prompter Padh in the 
test strain D-K2 (Figure 6). These observations indicate that 
alterations in gene expression in response to nutrient availability 
can be achieved through anti-regulation of expression between 
partners of a convergent gene pair, mediated by their 3'-UTR. 

We further tested the role of anti-regulation between convergent 
genes using a growth assay of the wild type strain (YL 1 C^^) and 
seven strains (D-Kl to D-K7, Table 2) genetically modified for 
the convergent gene pair ADEl and KIMS. In an adenine rich (SC) 
medium, the growth phenotype of all the strains largely reflects the 
change in the inoculation concentration (Figure 7, right panel). In 



the adenine-barren medium (SC-A), growth phenotype of the wild 
type strain is comparable to that in the rich medium, given that the 
ADEl gene is expressed at normal levels, as shown previously [22]. 
Growth of the strain with the KIN3 gene knocked out (D-K3) is 
comparable to that of the wild type strain, agreeing with the fact that 
expression level of the ADEl gene is comparable between these two 
strains (group III, Figure lA). However, growth is clearly repressed 
in the SG-A medium for strains in which ADEl expression was 
suppressed, either directly (D-Kl) or indirecdy by over-expression of 
its convergent partner gene KIN3, in situ (D-K2 and D-K4) or 
ectopically (D-K5). Ectopic over-expression of the KIN3 3' -UTR 
alone (D-K7) clearly repressed growth in comparison to the wild 
type, but no repression was observed with ectopic over-expression of 
the KIN3 ORF alone (D-K6). Combining this result with our earlier 
observations suggests that ectopic over-expression of eitlier the KIN3 
3' -UTR or of the entire gene suppressed the expression of its 
convergent partner gene ADEl, which in turn, led to suppressed 
growth in the nutrient limited medium. We conclude that the 
convergent organization of genes with overlapping 3'-UTRs in the 
yeast genome constitutes an eflFective mechanism for regulating gene 
expression and ultimately for controlling cell growth. 

Discussion 

The mechanisms of transcriptional regulation are highly con- 
served in eukaryotic species [23]. However, the yeast genome is 
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Figure 2. Dual-luciferase reporter assay of protein levels for two convergent gene pairs. Dual-luciferase reporter assay for the expression 
of the upstream genes of two pairs of convergent genes with overlapping 3'-UTRs. The reporter protein replaced the ORF of each upstream gene (A 
APTl or B ADEl), while the downstream genes were modified: (I) wild type control with no modification; (II) over-expression of the entire gene; (III) 
ectopic over-expression of the ORF alone; and (IV) ectopic over-expression of the 3'-UTR alone. Mean and standard deviation based on three 
replicates are shown for mRNA expression level (upper panel) and protein abundance (lower panel), measured in relative units compared to the wild 
type control, which was assigned a value of 1.0. 
doi:1 0.1 371/journal.pgen.1 004021. g002 
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Figure 3. Loss of anti-regulation in the absence of overlapping 3'-UTRs. Expression levels of the upstream and downstream genes of two 
convergent gene pairs with their overlapping 3'-UTRs removed. Mean and standard deviation based on three replicates are given for the expression 
levels: (I) for a baseline control, assigned a value of 1.0; (II) with inhibition of the upstream gene; or (III) with over-expression of the downstream gene. 
doi:1 0.1 371 /journal.pgen.1 004021. gOOB 



distinct in several ways [24]. First, by the complete absence of 
RNA-dependent regulatory systems (including miRNAs and 
siRNAs) present in higher eukaryotes. Second, by its highly 



compact structure, for example, over 20% of ORFs in the yeast 
genome are arranged in cir as sense and antisense gene pairs (S/ 
AS) that overlap and are transcribed from opposing DNA strands. 




Figure 4. Expression patterns of four convergent gene pairs during the yeast cell cycle. Expression of four pairs of yeast convergent genes 
with overlapping 3'-UTRs in cells from the wild type strain YL1 A, the a-mating strain of YL1C {Pi ADE1/KIN3, B AXL2/REV7, C APTl/UNGl and DSHMl/ 
YPT10). Gene expression was assessed 10 minutes from the start of the cell-cycle and at 20 minutes intervals thereafter. 
doi:1 0.1 371/journal.pgen.1 004021 .g004 
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Figure 5. Cell-cycle expression of the gene pair ADE1 and KINS, in convergent and tandem orientation. (A) Three strains tested: (I) wild 
type, in which ADE1 and K\N3 were in convergent orientation; (II) tandem URA3, in which KIN3 was replaced with URA3, and ADEl and URA3 were in 
tandem orientation; and (III) tandem KIN3, in which ADE1 and K\N3 were in tandem orientation. (B) Cell cycle expression of the three genes ADE1, KIN3 
and URA3 in the wild type strain. (C) Cell cycle expression pattern of the genes ADEl and KIN3 in tandem orientation. (D) Expression of the genes 
ADEl and URA3 in the tandem URA3 strain. 
doi:1 0.1 371 /journal.pgen.1 004021 .g005 



To explore the expression pattern of genes arranged in this way, 
we analyzed a complete set of over 600 convergent gene pairs with 
overlapping 3'-UTRs. The convergent genes showed a significant 
and consistent pattern of negative correlation in expression across 
a broad range of growth conditions, though this negative 
correlation could only explain a limited fraction of total gene 
expression variability. The mechanism underlying this widespread 
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in SC in SC-A in SC in SC-A 
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Figure 6. Anti-regulation of a convergent pair in different 
nutritional environments. Expression of the convergent gene pair 
{ADEl and KIN3) with overlapping 3'-UTRs in two yeast strains cultured 
in complete (SC) or adenine barren (SC-A) medium. Strain D-K1 has 
ADEl inhibited, while strain D-K2 has KIN3 over-expressed. 
doi:1 0.1 371/journal.pgen.1 004021 .g006 



transcriptional interference was explored in detail for four 
representative pairs of convergent genes. 

We analyzed the expression patterns of the four convergent 
gene pairs in wild type yeast and its derived strains with various 
genetic modifications. The data revealed that convergent genes 
regulate each other's expression such that an increase in expres- 
sion of either gene produces a decrease in expression of the other. 




SC-A SC-A 



Figure 7. Phenotypic effect of anti-regulation between a 
convergent pair in a nutrient-limited environment. Growth 
phenotype of the wild type strain (YLIC™^) and seven strains (D-K1 to 
D-K7; see Table 2) genetically modified for the convergent gene pair 
ADEl and KIN3. Strains were cultured in complete (SC) or adenine 
barren (SC-A) medium. 
doi:1 0.1 371/journal.pgen.1 004021 .g007 
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We showed that this pattern of 'anti-regulation,' or transcriptional 
interference occurs in single cells and does not require interaction 
between the proteins encoded by the convergent genes, since the 
effects persisted when the downstream gene was mutated so that a 
fully functional protein could not be made. Furthermore, we 
demonstrate that transcriptional interference between convergent 
genes is reflected in the abundance of the proteins encoded by the 
partner genes, showing that it has far-reaching effects in yeast cel- 
lular networks beyond the immediate changes at the transcript level. 

The functional importance of 'anti-regulation' in gene expres- 
sion was demonstrated using a time course analysis of the cell cycle 
for three pairs involved in this process. In general, a rise in the 
expression of one gene occurred in parallel with falling levels of 
expression of the partner gene. This led to peaks in expression of 
one gene corresponding with troughs in the expression of the other 
gene, at various stages in the cycle. Moreover, we demonstrate that 
this cell-cycle expression pattern remains even when the conver- 
gent orientation of the pair is artificially converted into a tandem 
orientation. We have noted that the Proudfoot group conducted a 
systematic survey of the regulatory mechanism of convergent 
genes in Schizosaccharomyces pomhe fission yeast. They demonstrated 
induced heterochromatinization for convergent genes occurring in 
a short period of Gl-S phase, leading to down-regulation of the 
genes. This heterochromatinization event is mediated by the RNA 
inference (RNAi) pathway. In addition, most RNAi genes are 
themselves in convergent arrangement, resulting in auto-regula- 
tion between the convergent pairs [12,25]. They modified 
arrangement of a pair of convergx'nt g(mes by inserting an extra 
gene between the com-ergent pair, causing loss of the Gl-S down 
regulation [25] . In contrast, here we have shown that modification 
of the convergent pair into a tandem arrangement does not result 
in loss of anti-regulation in budding yeast. This is likely due to 
fundamental differences in the biology of budding yeast (studied 
here) compared with fission yeasts [12,25], as well as to differences 
in experimental procedures. First, as in metazoa. Dicer, Argonaut 
and RdRP, which are essential in the RNAi pathway, are 
conserved in S. pombe fission yeast, but lost in S. cerevisiae budding 
yeast [26] . Additionally, RNAi mediated regulation of convergent 
genes in fission yeast occurs through heterochromatinization, 
which requires several jjrotcins including Clr4, RdRP, Tas3, Chpl 
and Swi6, for which there are no known orthologs in S. cerevisiae 
budding yeast [26,27]. These findings suggest there may be 
different mechanisms regulating the expression of convergent 
genes in S. cerevisiae budding yeast and S. pombe fission yeast. 
Second, we created a truly tandem arrangement of the gene pair 
by switching the orientation of one gene in the pair. Moreover, in 
the present study, the convergent gi-iu^s had intact 3'-UTRs even 
after their orientation was converted into tandem. Third, the 
present study has focused on only those convergent genes with 
overlapping 3'-UTRs. 

We have also demonstrated that the well-known environmental 
r(;spoiisc of yeast cells to changes in nutrient availahilit\" is mediated 
through changes in expression of both the environment-responsive 
gene ADEl, and its convergent partner, KIN3. We conclude that the 
convergent organization of genes in the highly compact genome of 
S. cerevisiae is functionally significant and speculate that it represents 
a mode of gene expression regulation that may compensate for the 
absence of RNA-dependent regulatory systems. 

The most widely accepted explanation for transcriptional 
interference between convergent gene pairs with opposite tran- 
scriptional direction is the transcriptional collision model [1 3] . The 
model proposes that during transcription, RNA polymerase progresses 
towards the 3' end of each gene, and so the two processes on 
opposite strands will eventually clash; this will result in a negative 



correlation in transcript abundance between the convergent genes 
[13]. One prediction given by the model is that the transcriptional 
interference between convergent genes should be released when the 
genes are no longer in their native convergent arrangement. For the 
four pairs of convergent genes studied here, we have demonstrated 
for the first time that transcriptional interference between the partner 
genes can persist even when one of those genes is no longer expressed 
from its original location. Specifically, we have shown that ectopically 
over-expressing either the entire downstream gene, or its 3' -UTR 
alone, leads to suppressed expression of the convergent partner gene; 
this suppression is dependent on the 3' -UTR and is not observed 
when the ORF alone is ectopically expressed. This indicates that 
the 3' -UTR sequences are both necessary and sufficient to mediate 
transcriptional interference of convergent genes in the S. cerevisiae 
genome. Furthermore, transcriptional interference, at least for 
these four gene pairs, cannot be explained by the transcriptional 
collision model and must involve other, unexplored mechanisms. 

Materials and Methods 

Definition of ORFs with overlapping UTRs in the S. 
cerevisiae genome 

The combination of high-resolution yeast transcriptome sequenc- 
ing data [17] and genome sequence data enabled us to predict both 
5' and 3' untranslated regions (UTRs) on a genome-wide basis. 
From all ORFs with either 5' (4,835) or 3' (5,212) UTRs, we iden- 
tified UTRs overlapping by at least 1 base pair. From the mRNA 
sequencing dataset, we identified 645 'convergent' ORF pairs with 
overlapping 3'-UTRs, 53 'divergent' pairs with overlapping 5'- 
UTRs and 65 'consistent' pairs with the same transcriptional 
direction. We also analyzed nascent transcript (NET-Seq) sequenc- 
ing data [18] and strand specific RNA (ssRNA-Seq) sequencing data 
[19]. These three datasets were summarized in Table S2. 

Yeast strains and growth media 

Haploid yeast strain YLIC was used as the wild type strain and 
is described in detail elsewhere [20]. Genetically modified strains 
shown in Table 2 were constructed from the wild type following 
PCR-based protocols as follows. The upstream inhibited (group I) 
or downstream over-expressed (group II) strains were constructed 
by inserting an inhibitor or promoter segment accordingly into the 
upstream region of the target ORF. The terminator Ttef and 
constitutive promoter Padh were amplified from the widely used 
plasmids pAG36 and pAG32 respectively. The downstream gene 
knock out (group III) or three ectopically expressed groups (V, VI 
and VII) were constructed through exchanging the relevant gene 
segments. The knocked out segments were extracted from the 
plasmid pAG36, and the downstream ORF and its 3' -UTR were 
amplified from the YLIC genomic DNA. For frame shift 
mutagenesis (group TV), one or two nucleotide bases were inserted 
into the coding sequence to generate premature termination in 
translation. In the terminator changed groups (VIII, IX and X), 
the 3 '-overlapping regions of the convergent pairs were replaced 
by a Tjvdh-Taoxi segment constructed through PCR-based 
fusion assembly [28]; the upstream gene inhibited (group IX) or 
downstream gene over-expressed (group X) strains were then 
constructed using the same methods used to construct group I and 
II strains. Finally, strains for the dual-luciferase assay (Table S4) 
were created as follows. Firstly, we used the reporter gene Flue 
(firefly luciferase) vectored in the plasmid pGL3 (Promega) to 
replace the non-overlapping ORF regions of upstream genes of the 
convergent pairs, and then inserted Rluc [Renilla luciferase) 
vectored in the plasmid pRL-SV40 (Promega) into the genome 
as the internal control. Flue and Rluc sequences can be found else- 
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where [29], together with the promoter and terminator sequences 
used for constructing the internal control. All the genetic modifi- 
cations constructed here were confirmed by Sanger sequencing. 
Primer sequences for molecular cloning are shown in Table S6. 

Unless specified, both wild t^pe and engineered strains were 
grown in the standard rich medium (YPD: 1% yeast extract, 2% 
polypepton, 2% glucose, plus 2% agar if necessary). In the adenine 
starvation environment, strains were grown in synthetic dropout 
medium (the adenine barren medium, SC-A: 0.67% yeast nitrogen 
base w/o amino acids, 2% glucose, 0.2% yeast synthetic dropout 
mixture without adenine, plus 2% agar if necessary). The synthetic 
complete medium (SC: 0.67% yeast nitrogen base w/o amino 
acids, 2% glucose, 0.2% yeast synthetic complete mixture, 2% 
agar if necessary) was used for the adenine rich environment 
control. Synthetic dextrose minimal medium (SD: 0.67% yeast 
nitrogen base w/o amino acids, 2% glucose) was used for the cell 
synchronization experiment. 

RNA extraction and real time quantitative PCR 

Total RNA was extracted according to the hot acid phenol 
method as described in [30], followed by DNase I (Promega) 
cleanup to remove contaminating genomic DNA as described in 
[31]. The fractionated RNAs were used in 3'-RACE and real-time 
quantitative PCR. After reverse transcription, 1 ^1 cDNA tem- 
plates were used for the quantitative PCR assay to compare expres- 
sion levels of relevant genes [32]. For every tested strain, we took 3 
independent clones as biological replicates for the PCR analysis. For 
each of the biological rephcates, there were 3 technical replicates. 
Expression level was presented as the ratio of normalized target 
concentrations (AACt), as suggested elsewhere [33,34]. 

3'-RACE sequencing 

Protocols for 3'-RACE were implemented as described 
previously in [35]. In detail, the fractionated RNA was reverse 
transcribed using the Oligo (dT) anchor primer. The cDNA was 
then amplified with the 3'-end PCR anchor primer and gene- 
specific primers. After gel-purification, a nested PCR with the 
second anchor primer was conducted as necessary. We inferred 
the sequence reads to be the 3'-end of transcripts whenever poly-A 
appeared. The length of 3'-UTRs was counted up to and includ- 
ing the last nucleotide base before the poly-A. The sequence of the 
anchor primers is available from the commercial protocol (3'- 
RACE System for Rapid Amplification of cDNA Ends, Invitrogen). 

Nascent RNA expression analysis 

Before newly transcribed RNA undergoes maturation process- 
ing (i.e. nascent RNA), the transcripts contain a segment which will 
be cleaved during the process of RNA maturation. Based on this 
structural feature of nascent RNA, together with the mRNA 
sequence and genomic sequence data, we designed nascent RNA 
specific reverse transcriptional primers to profile the expression of 
the nascent RNA. We tested for nascent expression of two pairs of 
com ergent genes in the yeast genome using their nascent RNA 
specific reverse transcriptional primers as follows: ADEl-nascent- 
RTl (5'-CACTGGCAAACAAGATATCG-3'), APTl-nascent-RTl 
(5'-ATATTACTAT TGCATATGCAGGTC-3'), KIN3-nascent- 
RTl (5'-AGAGACTGGCTTACTGCTAATAAG-3'), and UNGl- 
nascent-RT 1 (5 '-AAATGATATGTTTCACGTCCTG-3 '). 

Dual luciferase assay 

Luciferase assays followed the protocol described previously 
[36,37] using the dual luciferase reporter (DLR) kits (Promega). In 
detail, the tested cells were grown in rich medium (YPD, 30°C) 



until the logarithmic phase (ODeoo — 0.7-0.9). After washing and 
re-suspension in 1 xPBS, the cultured cells were maintained in 
100 |xl 1 X passive lysis buffer for 15 s. An aliquot of 5 |jl was then 
extracted from the buffer for scoring luminescence measurements 
with 25 [tl LAR II reagent from a Lumat LB 9507 (Berthold 
Technologies) set with 2 s delay time and 10 s measurement time. 
The same procedure was implemented after 25 |il Stop & Glo 
reagent was added. The protein expression level was recorded as 
the ratio of the firefly luciferase activity to the Renilla luciferase 
activity (Flue/ Rluc). For each test strain, at least three independent 
cultures were assayed. 

Cell-cycle synchronization 

The protocol for cell-cycle synchronization was implemented as 
described previously [38]. The a pheromone-responsive strain 
YLIA (MATa, harlA) was used to achieve cell-cycle synchroni- 
zation in exactly the same genetic background as YLIC, which 
differs only in the mating type. The BARl gene of YLIA was 
knocked out so that the strain would respond to a low density of a- 
factor [39,40]. The logarithmic phase cells (ODgoo = 0.5) were 
arrested and incubated for 1.5 hours once the ot pheromone level 
reached 50 ng/ml. The cells were subsequently released from 
arrest by pelleting after repeated washing with pre-warmed 
ddH20 (30°C). The cells were then suspended in pre-warmed 
SD medium with Pronase E (0.1 mg/ml Pronase E, pH 6.4, 30°C) 
for 10 min. Subsequently, 15 ml of the prepared cell samples was 
taken every 20 minutes over the next two hours (approximating a 
complete cell-cycle) for extracting RNA, while the cell mass was kept 
at 30°C. In total, expression levels were measured at 6 time points. 

Nutrient-dependent growth test 

Cells from the tested strains (D-Kl through to D-K7) were first 
cultured in rich medium (YPD) overnight, and then diluted to 
3x10'^ cells/ ml (OD^oo = 0.01). The diluted cultures were dropped 
on the adenine barren medium (synthetic dropout medium, SC-A) 
or the adenine rich medium (synthetic complete, SC) respectively. 
The cultures were further diluted in a gradient and 5 ^ll of every 
diluted culture was dropped on the test plates, which were 
incubated at 30°C for 48 hours. 

Profiling gene expression in single cells 

5 ml yeast cells of wild-type YLIC were grown to the 
logarithmic phase in SD medium. Cells were then washed twice 
with sterilized water, and re-suspended in 5 ml ultrapure water. 
Cells were counted using a hemocytometer and diluted with 
reverse transcription buffer (Invitrogen) to only one cell per 7 |il. 
7 |j.l ahquots were deposited into each well of a 384-well cell plate 
(Corning). We identified wells containing only a single cell (single- 
cell well) using an inverse microscope. The method used to profile 
gene expression in a single cell was slightly modified from the 
documented protocol [41], briefly described as follows. 1 |.il 
mixture of lyticase (2 Hg/jXl, Sigma), DNase I (1 unit/|xl, Promega) 
and RNase OUT (2 uint/ [d, Invitrogen) were added to the single- 
cell well, and these wells were incubated for 15 min at 30°C, then 
for 15 min at 37°C to lyse the cells. Another 1 ^ll mixture of 
Proteinase K (0.1 |J.g/|J.l) and human total RNA spike-in control 
(10 ng/|xl) was added to the single-cell lysatc- and incubated for 
10 min at 65°C. Each microliter of the reverse transcriptase 
mixture (1 mM dNTP, 5 |J,M oligo-dT, 5 |xg/|xl BSA, and 20 
unit/|xl Superscript III) was used to initiate the aforementioned 
single-cell reverse transcription. To enhance the template level 
for quantitative PCR assay, we first performed a 15-cycle nested 
PCR for KIN3 and ADEl genes, with human ACTl as the spike- 
in control. The nested PCR primers in 5' to 3' direction were 
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GCCAC.4ACATACGTCGGTACA and AGGATTTTTTCAA- 
TGTTTGTCAGC for KINS, TCTTCACCCCATCGACCAA 
and CAGTAAGCCAGTCTCTTAAAAATTGC for ADEl, GC- 
ACAGAGCCTCGCCTTT and CGTGCTCGATGGGGTAC- 
TTC for ACTl. A 0.5 |J,1 aliquot of PGR product was used as the 
template for the next round of RT-PCR, using the following primers: 
GCCACAACATACGTCGGTACA and GGGAGTATGGTT- 
GGTCCATCA for KINS, TCTTCACCCCATCGACCAA and 
GGGCAGGAGAGATGTTTTCG ior ADEl, GCACAGAGCC- 
TCGCCTTT and GTTGTCGACGACGAGCG for ACTl. 

Microarray datasets and analysis 

Microarray expression datasets were downloaded from the 
NCBI Gene Expression Omnibus (GEO) database under the series 
accession numbers GSE19213, GSE13684 and GSE5185 (http:// 

www.ncbi.nlm.nih.gov/geo). These datasets were collected from 
an S. cerevisiae wild type strain BY4741 cultivated in seven different 
environments (Table S3). The tested cells were harvested at the 
early log (or exponential) phase (ODgoo — 0.3-0.4). Total RNA was 
extracted and processed according to the manufacturer's instruc- 
tions (www.afFymetrix.com). The Yeast Genome 2.0 (YEAST 2.0) 
microarray, which containc'd 5,744 probe sets interrogating all 
annotated ORFs in the S. cerevisiae genome, was employed to 
profile transcript abundance of the ORFs. Transcript abundance 
was extracted from raw hybridization signal intensities of each 
probe set using the Robust Multichip Average (RMA) method 
implemented in R [42] . The expression value for each ORE was 
log2-transformed prior to analysis. 

Supporting Information 

Figure SI Status of the 645 convergent gene pairs with 
overlapping 3'-UTRs as predicted from mRNA-Seq data, in 
both nascent RNA sequencing data (A) and strand-specific RNA 
sequencing data (B). 
(TIF) 

Figure S2 Correlation in expression between convergent gene 
pairs with overlapping 3' - UTRs across seven growth environments. 
Box plots showing the expression levels of 645 convergent gene pairs 
with overlapping 3'-UTRs across sc\-cn growth conditions described 
in Table S3. Above each box is the Pearson correlation coefficient 
(r) between partner genes and the corresponding P value. 
(TIF) 



Figure S3 Mapped reads of convergent gene pairs APTl /UNGl, 
ADEl/ KINS, SHMI/TPTIO, and AXL2 / REV7 trom nascent RNA 
sequencing data (a-d) or from strand-specific RNA sequencing 
datasets (e-h). 
(TIF) 

Figure S4 Correlated expression of the convergent gene pair, 
KINS and ADEl, measured in single cells. 

(TIF) 

Table SI ORFs with overlapping 3 '-UTRs identified in the 

budding yeast S. cerevisiae genome. 

(DOC) 

Table S2 Summary information of RNA-Seq, NET-Seq and 
ssRNA-Seq datasets. 

(DOC) 

Table S3 Seven environmental treatments used in yeast 
microarray analysis studies. 

(DOC) 

Table S4 Genetically modified strains used in the dual luciferase 
assays. 

(DOC) 

Table S5 Sequence of the convergent gene pair, KINS and 
ADEl, in wild type and two genetically modified strains. 
(DOC) 

Table S6 Primers using in RT-PCR, 3 '-RACE, and molecular 

cloning. 

(DOC) 
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